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ABSTRACT 

Megiddo  introduced  a  technique  for  using  a  parallel  algorithm  for  one  problem  to  construct  an 
efficient  serial  algorithm  for  a  second  problem.  We  give  a  general  method  that  trims  a  factor  of 
C>(logn)  time  (or  more)  when  this  tedmique  is  applied. 
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1.  Introductloa 

We  solve  a  somewhat  unusual  sorting  problem  optimally;  it  is  described  in  section  2.  the 
sorting  problem  was  motivated  by  and  derived  from  its  application.  The  application  is  an 
improvement  of  Megiddo's  ingenious  technique  [Ml]  which  uses  an  efficient  parallel  algorithm  for 
one  problem  to  produce  an  efficient  serial  algorithm  for  a  second  problem.  The  general  idea  of 
Megiddo's  technique  is  as  follows.  Suppose  a  certain  problem  A  is  soved  in  T  time  units  on  a  /*• 
processor  machine.  And  suppose  a  serial  algorithm  for  problem  A  running  in  time  U  could  be 
applied  to  designing  a  serial  algorithm  for  problem  B  with  a  running  time  of  0{CU),  where  essen- 
tially the  algorithm  for  problem  B  carries  out  each  step  of  the  algorithm  for  A ,  taking  time  C  per 
step.  By  using  the  parallel  algorithm  for  A  as  a  serial  algorithm  we  would  obtain  a  serial  algo- 
rithm for  B  running  in  time  0(CTP).  Megiddo  showed  that  in  many  cases  problem  B  can  actually 
be  solved  in  time  O(CTlogP),  if  /*  is  not  too  large.  CXir  contribution  is  to  reduce  this  time  to 
0(CT)  in  many  cases.  We  describe  Megiddo's  technique  and  the  improvement  in  section  3.  Using 
the  improvement,  we  trim  a  factor  of  O(logn),  or  more,  from  the  running  times  of  eight  algo- 
rithms, which  we  describe  in  section  4. 

2.  The  Sorting  Problem 

We  describe  the  sorting  problem  as  a  game.  There  are  two  players:  the  sorter  and  the 
adversary.  The  game  is  played  on  a  sorting  network  for  sorting  n  items;  the  network  has  width 
n/2  and  depth  f(n).  The  idea  of  the  game  is  to  carry  out  a  sort  on  the  network;  however  only  the 
adversary  is  allowed  to  resolve  comparisons  (i.e.  to  determine  which  of  two  imputs  is  the  larger), 
and  the  adversary  is  uncooperative.  More  specifically,  the  two  players  alternate  their  moves, 
which  are  as  follows. 

(1)  The  sorter  assigns  weights  to  all  the  comparators  for  which  both  inputs  are  known  (the  active 
comparators).  Let  the  sum  of  the  weights  assigned  be  W  (the  active  weight).  Let  C  be  the 
comparison  corresponding  to  comparator  C.    K  C  is  assigned  weight  h-  we  consider  C  to 


have  been  assigned  weight  w  also. 

(2)     The  adversary  is  obliged  to  resolve  sufficiently  many  of  the  weighted  comparisons  so  that  the 
sum  of  the  weights  of  the  resolved  comparisons  is  at  least  W/2. 

The  sorter  wins  when  the  sort  is  complete.  The  aim  of  the  sorter  is  to  win  quickly.  A  turn 
consists  of  one  move  of  each  player;  we  show  the  sorter  can  win  in  0(logn  +/(n))  turns.  By  'sort- 
ing on  the  network'  we  mean  that  the  sorter  must  resolve  exactly  those  comparisons  that  arise 
when  the  network  is  used,  and  must  follow  the  ordering  of  the  comparisons  created  by  the  net- 
work (that  is,  if  an  output  of  comparator  D  is  an  input  to  comparator  C,  then  the  comparison  at  D 
must  be  resolved  before  the  comparison  at  C  is  attempted). 

We  play  the  game  in  phases,  starting  with  phase  1.  A  phase  will  consist  of  5  turns.  During 
the  ith  phase  each  active  comparator  at  depth  j  in  the  network  is  given  weight  4'~J.  We  prove  the 
following  invariant. 

Lemma  1:  At  the  start  of  a  phase  the  active  weight  is  bounded  by  n/2. 

Proof:  At  the  start  of  the  first  phase  we  have  n/2  active  comparators  at  depth  1,  and  all  other 
comparators  are  inactive.  So  initially  the  invariant  holds.  We  note  that  in  going  from  the  ith 
phase  to  the  i+lst  phase  we  increase  the  weight  of  every  active  comparator  by  a  factor  of  4,  so 
for  the  invariant  to  hold  it  is  sufficient  to  show  that  at  the  end  of  the  ith  phase  the  active  weight  is 
bounded  by  n/S. 

We  show  that  in  one  turn  the  active  weight  is  reduced  by  at  least  one  quarter.  Consider  an 
active  comparator  C  of  weight  w,  and  suppose  the  corresponding  comparison  C  is  resolved.  Then 
C  ceases  to  be  active,  and  up  to  two  comparators,  each  of  weight  w/4,  may  become  active.  So  the 
resolution  of  C  reduces  the  active  weight  by  at  least  w/2.  Let  the  active  weight  be  W.  In  one 
txim,  we  are  guaranteed  that  the  comparisons  resolved  by  the  adversary  have  combined  weight  at 
least  W/2.  Thus,  in  one  turn,  the  active  weight  is  reduced  from  W  to  at  most  3W/4.  On  repeating 
5  times  we  reduce  the  active  weight  from  W  to  less  than  W/4.  Substituting  W  =  n/2  we  observe 
that  at  the  end  of  a  phase  (5  turns)  the  active  weight  is  bounded  by  n/S.  □ 


We  now  show  that  this  process  terminates  reasonably  quickly. 

Lemma  2:  For  k^j+]J2\ogn,  during  the  kth  phase  there  are  no  active  comparators  at  depth 7. 

Proof:  Suppose  there  was  such  a  comparator.  It  would  have  weight  4'"-'&4^'^'°8''  =  n.  But  the 
total  weight  of  the  comparators  is  bounded  by  n/2.   Contradiction,  d 

Coronary:  The  sorter  wins  after  0(f(n)  +  \ogn)  turns. 

Proof:  After /(n)  +  l/21og/j  phases  there  are  no  active  comparators  (by  lemma  2).  That  is,  all  the 
comparisons  are  resolved,  so  the  sorting  is  complete.  Each  phase  consists  of  5  steps.  Hence  the 
sorter  wins  after  0(f(n)  +  ]ogn)  turns,  o 

Remark:  Lemma  1  depends  on  the  fact  that  the  indegree  and  outdegree  of  each  comparator  is 
bounded  by  two.  Lemma  1  still  holds  when  we  replace  the  bound  of  two  by  any  constant  bound, 
if  we  change  the  weights  appropriately  and  change  the  number  of  turns  in  a  phase.  When  we  have 
arbitrary  in-  and  out-degrees  it  appears  unlikely  that  lemma  1  holds. 

3.  Improving  Megiddo's  techniqae 

In  general  terms,  Megiddo's  technique  provides  a  way  to  search  a  partially  ordered  space,  of 
polynomial  size,  without  actually  constructing  the  space  (the  polynomial  may  be  large).  Instead  an 
implicit  binary  search  of  the  space  is  made.  We  use  a  sorting  algorithm  to  guide  this  search 
(Megiddo's  technique  is  not  restricted  to  sorting  algorithms).  Typically,  a  sorted  order 
corresponds  to  a  region  of  the  space  being  searched.  So  resolving  a  comparison  in  the  sort 
corresponds  to  reducing  the  size  of  the  space  in  which  a  solution  is  known  to  lie.  As  might  be 
expected,  a  single  comparison  is  expensive;  but  surprisingly,  for  some  problems,  several  comparis- 
ons can  be  batched  relatively  cheaply.  This  leads  Megiddo  to  use  a  parallel  sorting  algorithm  to 
guide  the  search,  for  in  the  searching  algorithm  he  can  batch  the  comparisons  that  are  performed 
simultaneously  in  the  parallel  algorithm. 


More  precisely,  suppose  we  have  a  fast  parallel  algorithm  for  one  problem,  problem  A  say, 
which  is  used  to  construct  a  fast  serial  algorithm  for  a  second  problem,  problem  B  say.  Further, 
suppose  problem  A  is  sorting  and  problem  B  has  the  following  unusual  features. 

(1)  Problem  B  can  be  solved  by  'sorting'  but  each  'comparison'  is  expensive,  that  is  it  takes  time 
C{n)  rather  than  time  0(1).   Typically  C(/i)  is  0(/i)  or  C>(nlogn). 

(2)  If  we  consider  m  of  these  'comparisons',  Ci,...,C„,  we  can  order  the  comparisons, 
C^(ijS  •  •  •  <C„(„),  in  the  following  sense.  (Wc  think  of  a  comparison  C  as  being  the 
question  'is  Cj<C2?',  which  has  answer  either  'yes'  or  'no'.)  An  answer  of  'yes'  to  C^yj 
forces  C„(i),...,C„y_j)  to  have  answer  'yes',  while  an  answer  of  'no*  to  C^g-.  forces 
C„y+n,...,C^(„j  to  have  answer  'no'.  Also  we  can  determine  the  relative  order  of  two  com- 
parisons fairly  quickly,  typically  in  time  0(1). 

Megiddo  does  not  describe  his  technique  as  applying  to  such  problems  for  it  is  more  general,  and 
in  fact  so  is  our  improvement.  Nonetheless,  we  use  this  formulation  both  for  simplicity  and 
because  many  of  the  problems  we  consider  have  this  form.  Next  we  describe  how  Megiddo  solves 
p)roblem  B  and  then  we  explain  our  improvement. 

Suppose  we  were  to  use  a  standard  efficient  sorting  algorithm  (running  time  O(nlogn))  to 
solve  B.  Then  we  would  obtain  an  algorithm  with  running  time  0(nlognC(n)).  Megiddo's  idea  is 
to  use  a  parallel  sorting  algorithm,  using  P(n)  processors  and  parallel  time  T{n).  At  each  time 
step  of  the  parallel  algorithm  we  have  up  to  P(n)  comparisons  to  resolve.  Instead  of  evaluating 
them  one  by  one,  we  solve  the  median  one.  This  immediately  resolves  half  the  comparisons.  We 
repeat  log(/*(n))  times  (-OQogn)  typically)  and  we  thereby  resolve  all  P(n)  comparisons.  So  we 
achieve  a  running  time  of  0(r(n)C(n)logn)  plus  overheads  for  nmning  the  parallel  computation 
and  finding  medians  of  sets  of  comparisons.  The  overheads  for  running  the  parallel  algorithm  are 
0{P{n)T{n)).  In  fact,  as  Megiddo  observed,  we  can  use  a  parallel  algorithm  which  runs  in  time 
T(n)  in  Valiant's  model  [V]  so  long  as  the  overheads  can  be  performed  efficiently  in  a  serial 
simulation.    To  find  the  median  comparison  we  use  the  fast  median  algorithm  [AHU,  pp. 97-99], 


running  in  time  0(P{n)),  assuming  ordering  two  comparisons  takes  time  (9(1).  Since  each  time 
we  halve  the  size  of  the  set  of  unresolved  comparisons,  the  time  taken  to  find  all  \og(P{n))  medi- 
ans that  we  need  is  also  0{P(n)).  So  over  T(n)  parallel  time  steps  we  take  time  0(P(n)T(n))  to 
find  medians. 

Thus  the  total  running  time  of  the  algorithm  for  problem  B  is 
0(/'(rt)r(n)  +  r(n)C(n)logn).  Typically,  P(n)  =  0(n\ogn),  T(n)  =  0(]ogn)  [P],  or  P(n)  =  0(n), 
r(n)  =  (logn)  [AKS].  (In  the  latter  case  the  constant  is  very  large.)  When  C(/i)  =  0(n)  (O(nlogn) 
respec.)  each  of  these  sorting  algorithms  gives  an  algorithm  for  problem  B  running  in  time 
0(n\og^n)  (O(nlog'n)  respec.). 

Our  improvement  is  to  trim  a  factor  of  logn  from  these  running  times.  We  use  the  network 
of  [AKS];  instead  of  performing  the  comparisons  as  described  above,  we  play  the  game  described 
in  section  2.  We  have  to  describe  how  to  simulate  the  adversary:  when  the  adversary  is  required 
to  resolve  a  weighted  half  of  the  comparisons,  we  resolve  the  weighted  median  comparison,  which 
by  (2)  above  immediately  resolves  a  weighted  half  of  the  comparisons.  (Finding  a  weighted 
median  of  n  items  takes  0(n)  time  [R].)  It  is  easy  to  see  the  time  taken  to  play  the  game,  apart 
from  the  comparisons,  is  O(nlogn).  (The  [AKS]  network  can  be  built  in  time  O(nlogn)  --  the 
constant  is  rather  large,  however.)  And  since  the  depth  of  the  [AKS]  network  is  OQogn)  we  need 
perform  only  O(logn)  comparisons.  So  for  C{n)-0(n)  we  have  a  running  time  of  O(nlogn),  and 
for  C{n)  =  0(n\ogn)  a  running  time  of  O(nlog'n). 

In  several  of  the  applications  problem  A  is  not  sorting;  instead  it  it  the  problem  of  finding 
the  minimum  or  of  finding  the  median.  Here  Megiddo  would  use  algorithms  having  O(loglogrt) 
[SV]  and  C>((loglogn)')  [CY]  parallel  steps,  respectively,  and  using  0(n/\og\ogn)  and  0{n)  proces- 
sors, respectively.  These  yield  running  times  of  (9(C(n)loglognlogn  +  n)  and 
0(C(n) (log!ogn)-^logn  +  ^(loglogn)-^),  respectively,  for  problem  B.  Instead,  by  using  a  sorting 
algorithm  we  achieve  a  nmning  time  of  0{C{n)]ogn  ■+  nlogn)  for  problem  B,  in  both  cases. 


4.  AppUcadoas 

We  discuss  ten  algorithms  to  which  we  have  applied  these  improvements.  They  are  drawn 
from  [M1,M2,M3,MT1,MT2,CSY].  To  give  the  reader  a  better  understanding  of  the  method  we 
describe  one  algorithm  in  detail:  fmding  the  weighted  1 -center  (problem  1,  below).  Our  descrip- 
tions of  the  other  algorithms  are  brief.  To  fully  understand  each  of  these  algorithms  the  reader 
should  refer  to  the  original  solution  (given  in  the  appropriate  reference). 

1)  Weighted  EUididean  1-Center  Problem  [M1,M2].  Megiddo  gives  two  solutions:  one  with  nm- 
ning  C>(n(logn)'(loglog/j)')  [M1,M2],  and  one  with  running  time  C>(nlog^n);  we  improve  the  latter 
solution,  obtaining  a  nmning  time  of  0(n\ogn).  Megiddo  and  Zemel  [MZ\  also  obtain  a  proba- 
bilistic algorithm  with  running  time  C>(nlogn).  In  this  problem  we  are  given  a  set  X  of  /»  points 
iaj,bi),  i=l,...,n,  in  the  plane,  each  a  with  positive  weight  w^,  and  we  seek  a  point  (jc,^)  so  as  to 
minimize  J(x,}')  =  max{wj[(x-ai)'^  +  (y-fc/)-^^'^:  i=l,...,n}. 

We  describe  the  solution  in  four  stages.  First,  a  method  to  determine  whether  a  given  point 
is  the  center.  Second,  a  method,  which  given  a  line  L,  determines  whether  the  center  is  on  L,  to 
the  left  of  L,  or  to  the  right  of  L.  Third,  a  method,  which  given  a  drcle  C,  determines  whether 
the  center  is  on  C,  inside  C  or  outside  C.  And  fourth,  a  method  to  find  the  center. 

To  determine  if  (x,y)  is  the  center  we  compute  the  point(s)  in  X  farthest  from  (x,y).  If  this 
is  a  unique  point  then  (x,y)  is  not  a  center  (for  by  moving  infinitesimally  towards  this  point  we 
reduce  the  maximum  distance  from  the  candidate  center).  If  there  are  several  farthest  points, 
(x,y)  is  a  center  if  and  only  if  there  is  no  straight  line  L  through  {x,y)  such  that  all  the  farthest 
points  lie  strictly  to  one  side  of  L .  If  there  is  such  a  line  L ,  then  the  center  is  on  the  same  side  of 
L  as  the  farthest  points.   This  test  can  be  carried  out  in  time  0{n). 

To  determine  on  which  side  of  a  line  L  the  center  lies  we  seek  the  point  p  on  L  minimizing 
d(p).  If  p  is  not  the  center  there  are  two  possibilities.  One,  the  farthest  point(s)  from  p  is  (are) 
on  the  same  side  of  L,  in  which  case  the  center  is  on  that  side  also.  Two,  the  farthest  points  from 
p  are  on  both  sides  of  L  (see  fig.l).    Since  p  is  the  point  on  L  minimizing  the  function  d(  ),  there 


are  farthest  points  to  both  sides  of  J ,  the  perpendicular  to  L  through  p.  Let  AT  be  a  line  through  p 
such  that  all  the  farthest  points  are  on  the  same  side  of  K  (there  is  such  a  line  since  p  is  not  the 
center).  Without  loss  of  generality,  assume  the  farthest  points  are  all  above  K  and  that  K  slopes 
downward  as  shown  in  fig.l.  Let  ^  be  a  farthest  point,  below  J .  The  center  is  no  farther  from  q 
than  is  p;  that  is,  the  center  lies  in  a  drclc  with  center  q,  radius  qp.  But  the  center  also  lies  above 
K\  we  deduce  the  center  is  to  the  right  of  L  (see  fig.l). 

So  we  seek  to  find  the  point(s)  in  X  farthest  from  p.  To  do  this  we  use  a  parallel  algorithm 
--  the  natural  tree  based  algorithm  running  in  time  C>(logn).  When  comparing  two  points  pj  and 
Pj  we  determine  a  drcle  (or  if  they  have  the  same  weight  a  straight  line)  such  that  p^  is  the  farther 
point  on  one  side  of  the  drcle  and  p,  is  the  farther  point  on  the  other  side.  The  circle  intersects  L 
in  at  most  two  points,  determining  three  segments.  To  determine  which  of  p^  and  p,  is  farther 
from  p  we  have  to  decide  on  which  segment  p  lies. 

Next,  we  describe  how  to  determine  whether  p  is  to  the  left  or  the  right  of  a  point  q  onL. 
We  determine  the  farthest  point(s)  at  ^.  If  q  is  not  the  center  let  AT  be  the  perpendicular  to  L 
through  ^;  if  the  farthest  point(s)  lie  strictly  to  one  side  of  K,  p  lies  on  this  side  also;  otherwise 
p=q.  This  test  takes  0(n)  time. 

Consider  the  first  n/2  comparisons  in  the  algorithm  for  finding  the  farthest  points  from  p; 
they  determine  at  most  n  points  on  L.  We  begin  a  binary  search  to  determine  between  which  two 
of  these  points  p  lies;  actually,  we  just  perform  the  first  two  steps  of  this  search.  So  we  either  find 
p  or  we  determine  that  p  lies  in  a  segment  not  holding  3n/4  of  these  points.  In  the  latter  case  we 
resolve  at  least  n/4  of  the  comparisons,  and  therefore  determine  that  at  least  n/4  points  are  not 
among  the  farthest  points  from  p.  This  takes  time  0(n).  We  now  seek  the  farthest  point(s)  from 
among  the  remaining  at  most  3n/4  points;  so  the  total  running  time  is  0(n). 

To  determine  on  which  side  of  a  circle  C  the  center  lies,  we  instead  show  how  to  solve  the 
more  general  problem  of  determining  on  which  side  of  a  closed  circulw  arc  C  the  center  lies.  A 
closed  circular  arc  is  the  perimeter  consisting  of  the  arc  and  the  striiight  line  joining  its  two 
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endpoints  (the  bounding  line).   We  maintain  the  following  invariant. 

Either  we  have  determined  on  which  side  of  the  dosed  circular  arc  C  the  center  lies,  or  we 
have  determined  a  subarc  C  of  C  such  that  the  center  lies  on  the  same  side  of  C  and  of  C . 

Initially  C  =  C.  We  aim  to  either  obtain  a  subarc  S  on  which  the  farthest  point(s)  is  (are)  unique 
(S  may  be  a  single  point),  or  to  find  on  which  side  of  the  dosed  circular  arc  the  center  lies.  We 
show  how  to  do  this  in  time  0(nlogn)  and,  if  the  first  case  holds,  how  to  determine  on  which  side 
of  S  the  center  lies. 

We  start  by  determining  on  which  side  of  the  bounding  line  for  C  the  center  lies.  If  it  is  on 
the  opposite  side  to  the  circular  arc  we  have  answered  the  query;  otherwise  we  continue.  As 
above,  we  seek  the  farthest  point(s)  in  X  from  some  point  on  S.  We  use  the  same  tree-based 
parallel  algorithm  as  above.  Each  comparison  of  two  points  p^  and  p,  determines  a  drde  Cp, 
such  that  pi  is  the  farther  point  on  one  side  of  C12,  and  P2  is  the  farther  point  on  the  other  side  of 
C12.  C12  intersects  C  in  at  most  two  points,  determining  at  most  three  subarcs  on  C ;  the  subarc 
to  which  we  restrict  our  attention  determines  the  result  of  the  comparison. 

Let  D  =  qrs  he  a  circular  arc  containing  S.  We  show  how  to  determine  either  on  which  side 
of  D  the  center  lies  or  wbidi  of  the  subarcs  qr  and  rs  contains  S.  Let  L^  and  L2  be  the  infinite 
lines  supporting  the  bounding  lines  for  qr  and  rs,  respectively.  We  determine  on  which  side  of  L^ 
and  Li  the  center  lies,  in  time  0{n).  There  are  four  possible  outcomes,  each  of  which  implies  one 
of  the  following  (see  fig. 2): 

(a)  the  center  is  outside  D 

(b)  the  center  is  inside  D 

(c)  qr  contains  5 

(d)  rs  contains  5. 

We  next  show  how  to  resolve  m/2  of  m  (<n/2)  comparisons  in  time  0(n).  The  m  comparis- 
ons determine  up  to  2m  points  on  the  arc  C .    Let  q  and  s  be  i-.e  endpcjiuts  of  C,  and  let  r,,  r-,, 


and  r^  be  the  quartiles  and  the  median  of  the  2m  points  (when  C  is  traversed  from  one  endpoint 
to  the  other).  By  applying  the  method  of  the  previous  paragraph  twice  we  either  determine  on 
which  side  of  C  the  center  lies  or  we  determine  which  subarc  among  qr^,  r-j-^,  r-fy^,  r-^s  contains  5. 
In  the  latter  cases  the  subarc  C"  containing  5  holds  only  mJl  of  the  2/n  points;  so  on  C"  (and 
hence  on  5)  the  outcome  of  at  least  m/2  of  the  comparisons  is  determined.  We  now  restrict  our 
attention  to  C" . 

By  using  a  game  similar  to  that  of  section  2,  in  O(logrj)  phases,  each  requiring  a  constant 
number  of  'comparisons',  we  either  determine  on  which  side  of  C  the  center  lies,  or  we  find  a 
subarc  5  of  C,  together  with  the  farthest  point(s)  from  5,  such  that  the  center  is  on  the  same  side 
of  S  and  of  C,  and  the  farthest  points  from  any  point  on  S  are  the  same.  In  the  latter  case  if  the 
farthest  point  is  at  the  center  of  C  we  deduce  the  center  is  inside  C;  otherwise,  let  p  be  the  point 
on  S  minimizing  the  function  (f( ).  It  can  be  found  in  a  farther  0{n)  time.  Let  i,  be  the  tangent 
to  S  at  /».  We  claim  the  center  cannot  lie  outside  5  and  between  the  tangent  and  the  semi-infinite 
extension  of  the  bounding  line  for  5,  which  touches  the  tangent  (see  fig.S).  Thus  by  determining 
on  which  side  of  L  the  center  lies,  we  determine  on  which  side  of  the  circular  arc  the  center  lies. 
The  overall  running  time  is  C>(/ilogn). 

Proof  of  Claim:  Suppose  the  center  lies  between  S  and  the  tangent.  Let  AT  be  the  line  joining  the 
center  to  p.  We  note  that  as  we  traverse  K  from  the  center  towards  p  the  value  of  the  function 
d(  )  is  strictly  increasing.  But  K  crosses  5  before  reaching  p ;  this  contradicts  p  being  the  minimum 
of  the  function  d()  on  5.  n 

Finally,  we  describe  how  to  find  the  center.  Again,  we  seek  the  farthest  points  from  the 
center.  We  use  the  same  tree-based  algorithm  as  above.  When  comparing  two  points  we  obtain  a 
drcle  C  (if  the  points  have  the  same  weight,  a  straight  line),  and  to  resolve  the  comparison  we 
have  to  determine  on  which  side  of  the  circle  the  center  lies. 

We  show  how  to  resolve  n/4  among  nfl  comparisons,  each  determining  a  circle.  Fu^st,  we 
obtain  a  horizontal  strip  in  which  the  center  lies  and  in  which  no  two  drcle  boundaries  intersect. 
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We  achieve  this  by  'sorting'  the  drcle  boundaries.  We  divide  each  drcle  into  its  left  and  right 
halves  (by  splitting  the  drcle  with  the  vertical  diameter).  Then  we  'sort'  the  n  semidrcles  using 
the  [AKS]  network.  How  do  we  'compare'  two  semicircles?  Let  Z-j  and  L;  be  the  horizontal  lines 
through  the  intersections  of  the  two  semicircles,  if  any.  A  comparison  reduces  to  determining  on 
which  side  of  Z,j  and  L^  the  center  lies.  So  a  comparison  takes  0(n)  time.  Using  Megiddo's  tech- 
nique and  our  improvement  we  'sort'  the  semidrcles  in  time  0(n\ogn).  When  we  have  found  the 
sorted  order  we  will  also  have  found  the  required  horizontal  strip. 

By  determining  for  two  left  semidrde  boundaries  on  which  side  the  center  lies  we  can  deter- 
mine to  which  side  (left/right)  of  3n/8  left  boundaries  the  center  lies  (out  of  n/2  left  boundaries); 
with  two  more  computations  we  can  determine  to  which  side  of  3n/8  right  boundaries  the  center 
lies.  Together,  these  determine  for  at  least  n/4  drdes  whether  the  center  is  inside  or  outside  of 
them.  Each  comparison  corresponding  to  such  a  cirde  is  thereby  resolved.  So  we  determine  that 
at  least  n/4  points  are  not  among  the  farthest  points,  in  time  O(nlogn).  We  continue  the  search 
among  the  remaining  at  most  3n/4  points.  Thus  we  find  the  farthest  points  in  time  O(nlogn). 
Given  the  farthest  points,  in  time  0(n)  we  determine  the  point  equidistant  from  them;  it  is  the 
center.  Thus,  we  find  the  center  in  time  O(nlogn). 

2)  The  Ham  Sandwkh  Theorem  in  2  EMmensions  [CSY].  In  the  discrete  ham  sandwich  theorem 
we  are  given  a  set  of  n  blue  points  and  a  set  of  n  red  points,  and  we  are  to  find  a  straight  line 
dividing  both  sets  into  equal  sized  halves.  The  algorithm  in  [CSY]  uses  a  parallel  algorithm  for 
finding  medians,  with  C(n)  =  0(n).  So  it  runs  in  time  0(nlogn(loglogn)-^.  Our  improvement 
yields  a  running  time  of  (9(«logn). 

3)  Finding  a  Point  In  the  Center  hi  2  Dimensions  [SCY].  A  center  of  a  set  of  n  points  has  the 
property  that  any  line  through  it  has  at  least  \n/3]  - 1  points  to  either  side.  [CSY]  first  solve  the 
problem  B,  of  determining  on  which  side  of  a  straight  line  the  center  lies,  in  time  0(n\og^n). 
Their  solution  uses  a  parallel  sorting  algorithm,  with  C(n)  =  0(n\ogn).  Our  method  improves  the 
0{n\og'n)  to  a  time  of  0(n\og^n).    [CSY]  next  solve  fl,.  ^  problem  of  'sorting'  n  lines;  the  solu- 
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tion  uses  a  parallel  sorting  algorithm  with  fl,  used  to  resolve  comparisons.  We  save  another  factor 
of  O(logn).  So  we  solve  fl,  ^  ^i™^  0(n\o^n)  rather  than  0(nlo^n).  The  main  algorithm  uses 
5;  O(logn)  times  so  our  version  of  the  algorithm  runs  in  time  0(/ilog"'n)  rather  than  0{n\o^n). 

4)  Finding  a  Point  in  the  Center  in  3  Dimensions  [CSY].  Here,  a  center  of  a  set  of  n  points  has 
the  property  that  any  plane  through  it  has  at  least  \n/4]  - 1  points  to  either  side.  The  algorithm  is 
similar  to  the  one  above.  [CSY]  first  solve  the  problem  B^  of  determining  on  which  side  of  a 
plane  the  center  lies;  their  solution  uses  n  parallel  sorting  algorithms,  running  in  parallel,  with 
C(n)  =  0(n^ogn).  Using  essentially  the  solution  for  problem  (2)  they  achieve  a  running  time  of 
0(n^o^n);  using  the  improvements  described  in  (2)  we  reduce  this  to  O(n'^og^n).  [CSY]  also 
solve  ^2,  a  problem  of  'sorting'  O(n')  planes;  we  save  a  factor  of  0(\ogn)  here,  also.  So  we 
reduce  the  running  time  from  0(n^og^°n)  to  0(n^og'n). 

5)  A  spanning  Tree  Problem  [Ml].  A  parallel  sorting  algorithm  is  used  to  sort  the  edges  of  the 
graph  with  each  comparison  costing  O(EloglogV).  So  Megiddo  obtains  a  running  time  of 
0(£(logV)'loglogV)  and  our  improvement  yields  a  running  time  of  (9(£logVloglogV). 

6)  A  Scheduling  Problem  [Ml].  A  parallel  sorting  algorithm  is  used,  with  each  comparison  cost- 
ing C>(nlogn).  So  Megiddo  obtains  a  running  time  of  0{nlo^n)  and  our  improvement  yields  a 
running  time  of  0(n]og~n). 

7)  The  Minimum  Ratio  Cyde  Problem  [Ml].  The  problem  is  to  find  a  cyde  in  a  network  with 
edge  costs  and  edge  times,  such  that  the  ratio  of  the  total  cost  to  the  total  time  of  the  edges  on  the 
cyde  is  minimized.  We  improve  Megiddo's  approach  No.2  [Ml,  p. 858].  He  uses  a  parallel  algo- 
rithm to  find  the  minimum  of  V  items,  each  comparison  taking  0(EV)  time.  This  algorithm  is 
iterated  0(]ogV)  times,  taking  time  0(V'logV  +  EV(logV)^og]ogV)  overall,  where  O(V'logV)  is 
the  time  for  computations  other  than  comparisons.  Instead  we  use  a  sorting  algorithm  to  find  the 
minimum;  together  with  our  improvement  of  Megiddo's  technique  this  reduces  the  running  time  to 
0(V'logV+£:v(logV)'). 
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8)  Max-Mln  Tree  k-Partldoning  Problem  [Ml].  The  problem  is:  given  a  tree  T  with  n  edges  and 
a  nonnegative  weight  associated  with  each  vertex,  delete  k  edges  of  7  so  as  to  maximize  the 
weight  of  the  lightest  of  the  resulting  subtrees  (formed  by  the  remaining  edges).  Megiddo  solves 
the  partitioning  problem  for  a  path  in  time  OinXo^n);  his  solution  requires  n  binary  searches  to  be 
done  in  parallel  on  a  set  of  n  items,  wliere  each  comparison  takes  time  0{n).  While  this  is  not  a 
sorting  algorithm,  we  can  still  apply  the  method  of  section  2  and  obtain  a  game  that  requires  us  to 
make  just  O(logn)  comparisons.  Using  this  game  we  solve  the  partitioning  problem  for  a  path  in 
time  0{n\ogn).  The  partitioning  problem  for  a  tree  is  solved  recursively,  using  the  solution  to  the 
path  partitioning  problem  to  put  the  pieces  together;  Megiddo  obtains  a  running  time  of 
O(nlog'n);  our  improvement  of  the  path  partitioning  algorithm  reduces  the  running  time  to 

9)  p-Center  Problems  [MTl].  See  [MTl]  for  the  problem  definition.  Their  solution  is  btiilt  on  a 
searching  problem  described  in  the  appendix  of  [MTl].  It  has  two  phases:  a  sorting  phase  and  a 
series  of  n  binary  searches  running  in  parallel.  For  both  of  these  the  comparisons  take  (?(n)  time. 
We  apply  the  method  of  section  2  to  the  sorting  phase  and  the  improvement  described  in  (8) 
above  to  the  binary  searches;  this  reduces  their  running  time  of  C>(nlog'n)  for  the  searching  prob- 
lem to  O(nlogn).  As  in  (8),  the  continuous  p-center  problem  is  solved  recursively,  using  the  solu- 
tion to  the  searching  problem  to  put  the  pieces  together.  Again,  Megiddo  obtains  a  running  time 
of  0{n\o^n)  for  finding  the  continuous  p-center;  our  improvement  to  the  algorithm  for  the  search- 
ing problem  reduces  this  to  (9(nlog-n). 

10)  Least  Distance  Lines  [MT2].  The  problem  is  as  follows.  Given  n  points  in  the  plane 
{x^,y^,...,{x„,y„)  together  with  positive  weights  w^,...,w„,  find  a  straight  line  L  so  as  to  minimize 

n 

2 >V/^(x, ,>-,;£,),  where  d  is  the  distance  function  from  L  relative  to  the  /j-metric.    Their  solution 

/-I 

uses  a  parallel  sorting  algorithm  to  sort  the  n  values  y^—ax^,  where  L  is  the  line  y  =  ax  +  b,  and 
both  a  and  b  are  unknown.  Thus  they  obtain  an  0(nlog'/j)  running  time.  Our  improvement 
reduces  this  to  0{n\ogn). 
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4.  Further  Remarks 

We  comment  now  on  the  efficiency  of  the  algorithms  described  above.  In  all  the  problems 
above  we  need  to  use  the  sorting  algorithm  based  on  the  [AKS]  network.  This  has  a  very  large 
constant  making  these  algorithms  impractical.  However,  we  could  implement  the  [AKS]  network 
probabilistically  and  then  the  constants  become  less  unreasonable,  while  the  order  of  the  expected 
running  time  is  unchanged.  In  addition,  when  a  median  is  required,  as  a  rule  an  approximation  to 
the  median  will  suffice,  so  if  we  are  using  a  probabilistic  algorithm  we  can  take  a  random  item 
instead.  This  may  not  improve  the  order  of  the  runtime,  but  in  practice  it  will  result  in  a  simpler 
and  more  efficient  algorithm. 

Unfortiinately,  we  cannot  apply  our  imprrovement  to  the  sorting  algorithms  of  Valiant  [V]  or 
Preparata  [P]  for  these  do  not  yield  sorting  networks  of  small  enough  depth.  An  efficient  parallel 
sorting  algorithm  {0(n)  processors,  O(logn)  running  time)  would  enable  all  the  algorithms 
described  above  to  be  implemented  more  efficiently. 

We  believe  that  other  seardiing  problems,  particularly  those  with  a  geometric  flavor,  will 
benefit  from  applying  these  techniques. 
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